function y = jacobi(x, fh)
    y = x;
    y(2:end-1) = (x(1:end-2) + x(3:end) + fh(2:end-1)) * 0.5;
end
